Dynamical preparation of EPR entanglement in two-well Bose-Einstein condensates 



B. Opanchuk 1 , Q. Y. He 1 - 2 , M. D. Reid 1 ^*, and P. D. Drummond 1 ^ 

1 Centre for Atom Optics and Ultrafast Spectroscopy, 
Swinburne University of Technology, Melbourne 3122, Australia 
2 State Key Laboratory of Mesoscopic Physics, School of Physics, 
Peking University, Beijing 100871, People's Republic of China and 
* pdrummondQswin. edu. au, * mdreid&swin. edu. au 

We propose to generate Einstein-Podolsky-Rosen (EPR) entanglement between groups of atoms in 
a two-well Bose-Einstein condensate using a dynamical process similar to that employed in quantum 
optics. The local nonlinear S'-wave scattering interaction has the effect of creating a spin squeezing 
at each well, while the tunneling, analogous to a beam splitter in optics, introduces an interference 
between these fields that results in an inter- well entanglement. We consider two internal modes at 
each well, so that the entanglement can be detected by measuring a reduction in the variances of the 
sums of local Schwinger spin observables. As is typical of continuous variable (CV) entanglement, the 
entanglement is predicted to increase with atom number, and becomes sufficiently strong at higher 
numbers of atoms that the EPR paradox and steering non-locality can be realized. The entanglement 
is predicted using an analytical approach and, for larger atom numbers, stochastic simulations based 
on truncated Wigner function. We find generally that strong tunnelling is favourable, and that 
entanglement persists and is even enhanced in the presence of realistic nonlinear losses. 



I. INTRODUCTION 

Entanglement between groups of atoms has been con- 
firmed experimentally, and recent experiments report the 
development of quantum correlated twin-atom beams [T]- 
[4|. This represents a first benchmark for investigations 
into multi-particle non-locality that could deepen our un- 
derstanding of the "classicality versus quantumness" for 
macroscopic objects [5j. So far, however, there has been 
no reported conclusive demonstration of stronger forms 
of quantum non-locality (such as violations of Bell in- 
equalities [6 , the Einstein-Podolsky-Rosen paradox [7H5] 
or steering [TUHH]) using mesoscopic groups of atoms, 
although there have been a number of theoretical pro- 
posals and studies of the correlation between spatially 
separated atoms [13 -20 . Highly efficient detection of an 
Einstein-Podolsky-Rosen (EPR) paradox for quadrature 
field amplitudes [5] El] > and loophole- free steering for pho- 
tons [12], has been been realized in optics. 

In this paper, our goal is to develop a strategy for gen- 
erating spatial entanglement between mesoscopic groups 
of atoms confined to the potential wells of an ultra-cold 
Bose Einstein condensate (BEC). Such atoms have been 
verified entangled |21j . but we seek to achieve an unam- 
biguous EPR paradox-steering type of entanglement, in 
which the entanglement between the atoms of different 
wells can be characterized, and readily extended to situ- 
ations involving the genuine multi-partite entanglement 
of groups of atoms at three or more sites [35]. Apart 
from the potential to test quantum mechanics in new 
regimes, this type of entanglement underpins many im- 
portant applications in the fields of quantum information 
and metrology [23H2S]- 

There are many possible strategies for the generation 
of such spatial multi-atom EPR entanglement. A re- 
cent experiment demonstrates EPR-type correlation near 
the coherent noise level using spin changing collisions 



[3J, and there has been a recent proposal to create spa- 
tial entanglement between two wells by direct adiabatic 
cooling to the ground state [57]. The most common 
method used in quantum optics however is to combine 
two squeezed single mode fields through a beam splitter 
(BS) El HH HH1 HI] ■ The method relies on a nonlinearity 
to produce squeezing in each mode locally, followed by 
a linear coupling transformation to create the entangle- 
ment between the two modes. 

Motivated by this, we explore in this paper a similar 
dynamical strategy applied to the BEC double poten- 
tial well system. Following Milburn et al [3U], we as- 
sume the atoms of each well can be modeled using a 
single mode approximation, and introduce the respec- 
tive boson operators, a and b. In this case, the S-wave 
scattering intra- well interactions, given by Hamiltonians 
H = hga^a 2 and H = Hgw 2 b , provide the nonlinear- 
ity at each well that generates a local spin squeezing, 
while the coupling or tunneling inter- we 11 term, mod- 
eled as H = Hnffib + aw), provides the linear beam- 
splitter transformation [3"Tl 13"2"] that generates inter-well 
two- mode entanglement. 

We consider in fact two dynamical strategies: in the 
first the local nonlinear and nonlocal tunnelling processes 
act sequentially; in the second they act simultaneously. 
While the second strategy is likely to be more practical, 
we analyze the sequential case first, in the earlier sections 
of the paper, because it allows a full quantum solution 
within the constraints of the two-mode model. 

We show that substantial entanglement can be gener- 
ated in both cases, provided parameters are optimized. 
An analysis of what is accessible experimentally indicates 
that this entanglement could be within reach of current 
set-ups, though the realization of EPR and steering para- 
doxes presents a greater challenge, and may require more 
sophisticated experimental procedures. 

Perhaps surprisingly, we find that a large amount of en- 
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tanglement requires using large numbers of atoms. This 
result is verified using both a full quantum analysis and 
the truncated Wigner function [33] which becomes valid 
as 1/iV 3 / 2 — » 0. The manifestation of a continuous vari- 
able (CV) EPR paradox has been shown for large optical 
amplitudes (34], but may have been thought impossible 
to realize for large groups of atoms. The prediction is 
consistent with those of Ferris et al [TS] and who study 
EPR entanglement in BEC four wave mixing, though in 
our case an atomic homodyne detection [3J is not used, 
the EPR observables being the local Schwinger spins that 
can be measured via Rabi rotation and atom counting. 

The paper is organized as follows. We first summarize 
in Section II the meaning of Einstein-Podolsky-Rosen en- 
tanglement, and outline how it is to be detected. The 
dynamical solutions for the entanglement via the two- 
step method are presented in Sections |III| and |IV| The 
truncated Wigner function simulations and the results 
for the inter-well entanglement via simultaneous evolu- 
tion are treated in Section [V] An analysis modeling cur- 
rent experimental regimes, including the effect of nonlin- 
ear losses, indicates that the EPR entanglement is robust 
against the expected decoherence effects. 



II. EPR ENTANGLEMENT 

The original CV EPR paradox (7] considers correla- 
tions between the positions and momenta of two parti- 
cles emitted from a source. With optical or atomic Bose 
fields, one can define the quadrature phase amplitudes of 



two spatially separated field modes, as Xa 



■a) /I, 



and P A = (fit - a)/2i, and X B = (&t + b)/2, and 
Pb = (b* — b)/2i. These are analogous to the position 
and momentum in the two-particle system. 

The EPR paradox arises when both Xa and Xb, and 
Pa and Pb, are maximally correlated, so that measure- 
ment of Xa enables exact prediction of the outcome for 
measurement of Xb , and the measurement of Pa enables 
exact prediction of the outcome of measurement of Pb- 
EPR argued that, assuming no "spooky action at dis- 
tance", the action of measuring Xa could not "create" the 
result for Xb- They then concluded, since the result for 
Xb can be predicted without disturbance of that system, 
that the result for the outcome is predetermined. Since 
the predetermination of both Xb and Pb is without un- 
certainty, there can be no equivalent local quantum state 
interpretation. Hence, EPR argued that quantum me- 
chanics was incomplete. The premises assumed by EPR 
are often termed "local realism", and the EPR paradox 
can be thought of as a demonstration of the incompat- 
ibility between "local realism" and the "completeness of 
quantum mechanics". 

EPR's argument applies when one observer (Alice) can 
make precise predictions for the outcome of measure- 
ments made by a second, distant observer (Bob). The 
key issue for the EPR paradox is that Alice can infer a 
result for either of two of Bob's conjugate observables, 



by measuring locally on her system. The EPR paradox 
arises when the accuracy of her inferences would violate 
quantum mechanics, if she could infer results for both 
conjugate observables simultaneously. This demonstra- 
tion of the EPR paradox is most simply achieved by 
comparing the conditional variances for Alice's measure- 
ments with the variances of the Heisenberg uncertainty 
relation [g]|35]. 



A. Entanglement criteria 

The original EPR paradox focused on states that 
showed correlation and anti-correlation for position and 
momentum respectively Duan et al and Simon [361 137| 
showed that entanglement between modes a and b is con- 
firmed if 



D = A 2 (X A - X B ) + A 2 {P A + P B ) < 1. 



(Ill) 



In the case ( II. 1 ) , the 1 arises from the commutation re- 
lation [a, fit] = 1, and reflects the quantum noise associ- 
ated with the four observables. A more general criterion 
is the product form: entanglement is confirmed when 
A(X A - X B )A(P A + P B ) < 1/2 El [35]. 

Quadrature phase amplitudes are measured using local 
oscillator methods, in which a strong field interferes with 
a signal field, using a beam splitter. A full analysis of 
a local oscillator measurement shows that it is actually 
equivalent to a Schwinger spin measurement, once the 
local oscillator is accounted for. It is useful to consider 
entanglement measures that have been developed for spin 
measurements. In particular, one can show entanglement 
using the spin version of ( II. 1 1 [29 , 



^{Ji t J#) + ^{Jl ±J Y B )< \(Ja)\ + \(Jb)\, (n.2) 

and also the Heisenberg-product entanglement crite- 
rion 138 



av*-j»).a2(j; 



S+it/2 
A 



J 



S+it/2 
B 



)< 



(Jl) 



(II.3) 



B. EPR paradox steering criteria 



While entanglement as confirmed by ( II. 1 



is necessary 

for the EPR paradox, it is not sufficient. To quantita- 
tively demonstrate the paradox, in the style constructed 
by EPR, the level of correlation in Alice's predictions 
must be compared with the quantum limit for a local 
state that might predetermine Bob's statistics. Thus, for 
the EPR paradox, the relevant quantum noise level is 
that of one observer, B, alone. 

The EPR pa radox confirms "steering", whereas the en- 
tanglement of (II.ll does not. Steering has been estab- 
lished as a distinct form of non-locality [TU]- The EPR 
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paradox and steering types of entanglement provide a dis- 
tinct resource for quantum information with applications 
not achievable for arbitrary entangled states |25| . 

An EPR paradox signature has been formulated in 
terms of conditional variances [3S]: thus EPR entangle- 
ment is observed when 



A{X B \X A )A(P B \P A ) < 1/4. 



(II.4) 



Since the precise choice of measurement at mode A is not 
important, only the inference, this criterion is sometimes 
more generally written as AinfA^Amf Pb < 1- 

More recently, it has been pointed out that the EPR 
signature is achieved once the entanglement variances be- 
come small enough |5J [TjJ]. EPR entanglement is de- 
tected, if either 

D = A 2 (X A - X B ) + A 2 {P A + P B ) < 1/2, (II.5) 

or A(X A - X B )A(P A + P B ) < 1/4. These criteria are 
special cases of (II. 4 1, and may not give the optimal mea- 



surement strategy for obtaining an EPR paradox, but are 
useful in many practical cases where only (II. 5 1 is mea- 
sured. 

Spin EPR measures can also be obtained, just as with 
entanglement per se, using the inferred Heisenberg uncer- 
tainty principle approach. From the Heisenberg spin un- 
certainty relation AJ B AJ B > \(J§)\/2, one can derive 
several simple spin-EPR criteria [9,29,39. For the case 
of a large non-zero mean spin, taken for definiteness in 
the z~ direction, an EPR paradox is demonstrated when 



A inf jf A inf j£ < ~|<Jf>|. 



(II.6) 



This criterion for the EPR paradox can be expressed in a 
particularly useful, though less general form, as follows: 



A2(4- 5 J£).A%7 



6+ir/2 
A 



9'J\ 



e+TT/2, ^ \(J B )\ 



< 



(II.7) 

In this version, the simplification has been made that the 
inference is given using a fixed linear approximation with 
constant gain g. 



III. DYNAMICAL PREPARATION OF BEC 
EPR ENTANGLEMENT 

We now turn to examining how to generate EPR en- 
tanglement in a BEC. A two- well BEC system can be 
modeled by the two-mode Hamiltonian [30 : 



H/h = k (a r b + aP 



tfataa + PPbb) . (III.l) 



Here k is the conversion rate between the two compo- 
nents, denoted by mode operators a and b, and g is the 
nonlinear self interaction coefficient proportional to the 
three-dimensional S-wave scattering length, 0,30- The 




Figure 1. Creation of spatial entanglement between modes 
and bi at different sites. The optical scheme is depicted in (a) 
and the BEC two- well proposal in (b). Local spin squeezing 
can be produced when the two modes at each well evolve for a 
time t under the nonlinear Hamiltonian ( |IV.1| | . Depending on 
the exact configuration, there can be coupling terms between 
the two modes at each site. The well height is controlled to 
allow a tunneling cross-interaction between the modes. This is 
equivalent to the modes interfering using a beam splitter, and 
the effect is to generate an inter- well spatial entanglement. 



first term proportional to n describes an exchange of par- 
ticles between the two wells (modes) in which total num- 
ber is conserved. This linear term is equivalent to that 
for an optical beam splitter |311 13"2"] . The Hamiltonian 
model applies to other two-mode bosonic systems includ- 
ing optical cavity modes or superconducting wave-guides 
with a nonlinear medium, as well as other types of two- 
mode atomic systems. However, we will use a two- well 
picture to illustrate the gedanken-experiment we have in 
mind, even though implementations may not be exactly 
in this form. The two-well model and system has been 
studied extensively in relation to macroscopic superpo- 
sition states |32l |40| and ultra-sensitive interferometric 
measurement |41l |4"2~] . 

If, in addition, we allow two internal spin components 
per well, the Hamiltonian becomes: 



H/h = ^ K i &f A + 9 X] h 



(III.2) 

Here we consider two internal modes at each EPR site A 
and B. with four modes in total, as shown schematically 
in Fig. [I] The local modes can be independent (in which 
case local cross couplings gij are zero), or not indepen- 
dent, as where the modes are coupled by the BEC self 
interaction term |42| . The coupling constant is propor- 
tional to the three-dimensional S-wave scattering length, 
so that gij cx ay , as in the two-mode case. 

We will describe our results using the equivalent di- 
mensionless Hamiltonian 



H = ^ iiialbi 



1^2 9ij a\a]ajCi t + 



with dimcnsionless coefficients k, 



{ai^bij, (III.3) 



Ki/gnN A , gtj 
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9ij I ' 9\\N A-, where Na is the initial total boson num- 
ber in well A, and a corresponding dimensionless time 
t = guNj^t. For definiteness, we will choose ratios 
of nonlinear couplings gij to correspond to known S- 
wave scattering lengths of 87 Rb, (between |1) = \F — 
1, 771 p = +1) and |2) = \F = 2, nip = —1)) at spe- 
cific magnetic field strengths near a Feshbach resonance 
B = 9.105 G. Namely, we use B = 9.086 G with corre- 
sponding ai2 = 107. 8ao, and B = 9.116 G with corre- 
sponding ai2 = 80.8ao (43], where ao is the Bohr radius. 
Intra-species scattering length are constant and equal to 
an = 100. 4ao and 022 = 95.5ao. Other results will be 
given as well. 

We propose to generate EPR states by preparing the 
system in a multi-mode coherent state and then allow- 
ing evolution according to the Hamiltonian (III.ll. Two 



strategies are compared. Firstly, where a manipulation of 
tunneling (e.g., by changing the potential barrier height) 
occurs, so that the nonlinear and linear terms are applied 
sequentially, and secondly, where the nonlinear and lin- 



ear tunneling terms act simultaneously, as in ( III. 1 III. 3 1 



The second strategy is likely to be more readily imple- 
mented, and is solved for in Section [V] via a truncated 
Wigner function method, with the inclusion of further 
effects such as nonlinear losses. 



In the two-step strategy, local squeezing is first gener- 
ated for each mode a and b via the nonlinear Hamiltonian 
in the absence of tunneling, and subsequently, a strong 
tunneling interaction provides a linear mixing that gen- 
erates an EPR entanglement between the two new modes 
c and d. The strategy is depicted in Fig. [T] The tech- 
nique is similar to that investigated experimentally in 
fibre-optics entanglement An important difference 

is that the fiber experiment used time-delayed pulses and 
dispersion to eliminate interactions, that is local cross 
couplings, between the components. This is not read- 
ily feasible in BEC experiments, although Feshbach res- 
onances can achieve this to a limited extent. In the next 
Section, we present full quantum solutions for this two- 
step strategy. 

We briefly remark on our choice of initial conditions. 
In experiments with multiple wells, it is generally possi- 
ble to prepare condensates with relative phase coherence, 
provided that tunneling is strong enough during the evap- 
orative cooling process. The overall phase is random, 
with a total number uncertainty that is typically at least 
Poissonian. This quantum state is therefore well repre- 
sented by a mixture of coherent states with a random 
over 11 phase. However, since none of our results depend 
on the overall phase, it is sufficient to consider just a sin- 
gle overall coherent state with an arbitrary phase. This 
is a low-temperature limit, which will develop additional 
fluctuations as temperatures are increased to the critical 
temperature. 



IV. STRATEGY I: TWO-STEP DYNAMICAL 
ENTANGLEMENT GENERATION 

A. Generation of local spin squeezing 

In the first step, squeezing is generated locally via a 
nonlinearity, given by the Hamiltonian 



i,3 



(IV.l) 



The initial state is a product coherent state for each 
mode: \a / 2) a Aa / y/2) ao . This models the relative co- 
herence between the wells obtained with a low inter- 
well potential barrier, together with an overall Poissonian 
number fluctuation that is typically found in an experi- 
mental BEC. We explain in the Appendix how to calcu- 



late the nonlinear quantum dynamical solutions of (IV.l I 



These are exact calculations, provided the original multi- 
mode Hamiltonian can be reduced to simple one or two- 
mode forms. 

Let ai,ci2 be operators for the two internal states at 
well A, and 61, 62 operators for two internal states at well 



B. Here Na = a\a,2 



a\a\ and Nb 



b\bi + b[bi are 



the atom number operators of these modes in each well. 
We introduce a phase-rotated Schwinger spin operator 
measurement that can be performed at each site. For 
site A, we define 



J a = 2 («2 a i e 



1 



Jl = 



2i 



a 2 aie lAe - a\a 2 e- lAe 



(IV.2) 



Ja = 2 (4^2 - a[ai 



where A9 = 9 2 — 0\ is the phase shift between the mode 
1 and mode 2. There is also an analogous definition at 
B with phase shift Atfi — 4>2 — <Pi- 

We select the phase shift to ensure ( J Y ) ^ 0, and the 
Schwinger spin operators orthogonal to J Y are given by 



j e = cos(d)J z + sin(6»)P 



(IV.3) 



all of which have the property (J ) — 0. Here Ad = 
A<j) = 7r/2 — a, and a is time dependent given by the 



character of (a 2 ai) 



This plane contains 



an infinite family of maximally conjugate Schwinger spin 
operators, generally given by J e and J e + 7T / 2 which obey 
the Heisenberg uncertainty relation 



A je A j(e+*/2) > 1(^)1/2. 



(IV.4) 



Quantum squeezing occurs when the variance in one con- 
jugate observable is reduced below the Heisenberg limit. 
Thus, 

A 2 J e < \(J Y )\/2 (IV.5) 
is said to be a "spin squeezed" state |45!l4"7] . 
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Figure 2. (Color online) Squeezing of local Schwinger spin op- 
erators versus interaction time r. The plot shows the squeez- 
ing A 2 j e /n (solid lines), and A 2 J e+w/2 /n (dashed lines), 
where the shot noise level is no = |(J y )j/2. Squeezing is ob- 
tained when 5" = A 2 j e ' e+n/2 /(\{J Y }\/2) < 1. Inset shows 
the optimal phase choice 6 for squeezing. The dimensionless 
coupling parameters correspond to 87 Rb atoms at magnetic 
field B — 9.116 G (blue lines), with corresponding scatter- 
ing lengths, as explained in the text. We also give results 
for the case without local cross couplings (red lines), i.e., 
512 = 0, 322 = Here N A = 200. 



Figure [2] shows the prediction for dynamical spin- 
squeezing, according to the solutions explained in the 
Appendix A, based on the Hamiltonian Eq. (IV. 1 ). Here 



we have assumed that a\ , b\ and ai , bi are initially in co- 
herent states. For simplicity, we suppose that the initial 
state is prepared in a four-mode coherent state by using 
a Rabi rotation: = I^WI^U^)^)^. After 
preparation, we assume that the inter- well potential is 
increased so that each well evolves independently, to give 
the solutions. 

We have considered the conditions required to obtain 
the best squeezing of Schwinger spin operators by opti- 
mizing the phase choice 9 (inset of Fig. [2]). We set 

dA 2 J e 



d 



89 



[cos 2 9A 2 J Z +sin 2 0A 2 J x 



09 1 

+2cos9sm9(J z , J x )] 
= 2cos(260(J z , J x ) - sin(26»)(A 2 J z - A 2 J X ) 
= (IV.6) 

and therefore obtain as the optimal squeezing angle 



1 



■ arctan . 
2 \ A 2 Jf - A 2 J? 



(IV.7) 



where 



(Ji,Ji) = 2 ( Wt) + (JiJX) - HJtKJi) 



The squeezing value in this case is: 

\{Ji)\n 

where 



A 2 T0,eW2 



(IV.9) 



A 2 J^ = cos 2 6» A 2 j| + sin 2 6> A 2 +2 sin 9 cos 9(JaJa) ■ 

(IV. 10) 



B. Producing the spatial entanglement 

The next step, after generating the local squeezing at 
each well, is to decrease the inter-well potential for a short 
time, so that it acts as a controllable, non-adiabatic beam 
splitter [35], to allow interference between the wells; other 
methods of obtaining an effective beam-splitter could also 
be feasible. 

Entanglement is generated by the interference of two 
squeezed states on a 50 : 50 beamsplitter with a rela- 
tive optical phase of <p (Fig.JTJ. The resulting entangled 
modes are labelled by c\^. and dip. Schwinger spin op- 
erators Jq/ d are denned for these modes, in accordance 



with Eq. ( |IV.2| . We note that f/ +7r/2 , f/ +7T/2 are 
measurable locally, in the style necessary for an EPR 
experiment, by using Rabi rotations and number mea- 
surements. 

The input-output relations for the Schwinger spin op- 
erators are given by c = id + re lip b, and d — tb — re~ lip a 
with the amplitude reflection and transmission coeffi- 
cients being denoted t and r. Here we use ip — ir/2 and 
r = t = 1/V2. Thus: 



a + ib ? b + id 
d = 



V2 



V2 



(IV. 11) 



The variances A 2 (j£ - J e D ) and A 2 (J* + " /2 + f+* /2 ) 



can both be small, so that 
S*_ = A 2 (J, 



'c 



Jb)< ^(\(Jc)\ + \(Jd)\) (IV.12) 



and 



5 + = A 2 (J^ /2 



J 



O+tt/2 
D 



)<^(\(Jc)\ + \(Jd)\), 



the degree of variance reduction being limited only by 
the amount of squeezing in the input modes. This is the 
signature of EPR entanglement, in accordance with the 
entanglement criteria (II. 2 1 and (II.3I. In fact, 



(IV.8) 



A 2 (JS T J e D ) = cos 2 6A 2 (J§ T J%) + sin 2 9A 2 (J X T J x ) 
+ cos 9 sin 6[(J§ qp -?E, Jc T J*) 
+ (J x T J x Jz T jZ)}. (IV.13) 

A similar expression can be given for A 2 ( J B C ± The 
solutions show a reduction in these variances, for suitable 
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(Color online) Inference squeezing: 
— Jd)/uo (solid lines) and S+ = 



\A 2 {J e c +lr/2 + J* +7r/2 )/n J (dashed lines) after the 
beamsplitter interaction, where the shot noise level is 
no = (\{Jc)\ + \(Jd}\)/2. Here N = N A = N B = 200. The 
parameters are the same as Fig. [2] The red lines show the 
result without local cross couplings, i.e., gi2 = 0, <?22 = <?ii- 



parameters, to indicate that Alice can infer Bob's observ- 
able Jq (by measuring j^,) with increasing accuracy, as 
A 2 ( J e c — Jjj) — > 0. Similarly, she can infer the conjugate 



observable J^, +7r ^ 2 (by measuring J^ ) ~ t " 7r/Z ) to an increas- 
ing accuracy, as A 2 ( J^ +7T ^ 2 + J^, +7r ^ 2 ) — > 0. We call these 
variances "inference variances" to remind us of their role 
in the EPR paradox. Ideally, we want to find a regime 
for which both inference variances become very small. 

It is convenient to express the output spin operators 
in terms of the inputs: we find 



fz _ fz 



J c + J n — 

tx 



D 

fX _ 
J n — 



■ t [a 2 b 2 — b\a 2 — a\b\ + b\ai], 
[a 2 a 2 — a\di + b\b 2 — b\bi], 
[e^-^aih-e^-^blar 



-e^-^a\b 2 + e-^-^b\a 2 ], 



J A + J R J, 



c 



J 



D > 



(IV. 14) 



Figure 4. Entanglement (-E pro duct < 1) based on the crite- 
rion in product form (IV. 15 1. EPR paradox entanglement is 
obtained when i? pr oduct < 0.5. The solid and dotted lines 
stand for N = N A = N B = 200 and for N = 2000 with cou- 
plings corresponding to B — 9.116 G; while the dash-dotted 
and dashed lines assume no local cross couplings, i.e., gi2 = 0, 
gn = ff22, for N = 200 and for N = 2000. 



criterion Eq. (II. 3 1: 



-^product — / * \ < 1. 

(l(^>l + K^)l)/2 

(IV.15) 

After using the beam splitter, entanglement can be de- 
tected as ^product < 1, as shown in Fig. [4] by the solid 
curve, which assumes the couplings between spins found 
at the rubidium Feshbach resonance. Note that, consis- 
tent with the results found in previous studies of entan- 
glement in the ground state, the dashed curve of Fig. [4] 
shows that no cross couplings, i.e., gi 2 = 0, gives much 
better results still. This would require spatially sepa- 
rated condensates for each spin orientation, in order to 
eliminate cross couplings, as recently demonstrated by 
using magnetic gradient techniques jJS]. Fig. [4] reveals 
improvement in the entanglement, as one increases the 
number of atoms in the condensate. 



C. EPR paradox entanglement 



where we have used #2 — 4> 2 = #i — 4>i = due to the 
symmetry of a and b, and 8 2 — 4>\ = 2 — Q\ = tt/2 — a as 
introduced for Eq. (IV.3I. The solutions for the inference 

and A 2 ' 



d+ir/2 
C 



J D 1 ) are 



squeezing A 2 (J£ - J 6 D ) and A 2 (J' L 
shown in Fig. [3] We note that unlike for most EPR states, 
the inference variances are asymmetrical. 

The spin orientation measured at each site can be se- 
lected independently to optimize the criterion for the 
state used. One can then show EPR entanglement via 
spin measurements by using the product entanglement 



The entanglement predicted for the dynamical scheme 
is strong enough that it reveals an EPR paradox (steer- 
ing) non-locality. This level of entanglement is reached 
when .Eproduct < 0.5 (as shown in Section UU, which oc- 
curs in the presence of local couplings for the larger atom 
numbers V = 2000, as shown in Fig. [4] The EPR para- 
dox entanglement can be obtained for lower atom num- 
bers N ~ 200 when local couplings are non-existent. 

Next , we will examine the predictions for the more sen- 



sitive EPR criterion Eq. (II.4I, which involves measure 



7 



ments of the conditional variances. The EPR argument 
is based on an accuracy of inference, that an observer at 
D can predict the result for an observer at C, to a 
certain measurable level of uncertainty. A simple way to 
determine an upper limit to this uncertainty is to use a 
linear estimate gJ^, based on the result J D for measure- 
ment at D. Then we arrive at the EPR paradox criterion 



of Eq. (11.71 



E 



EPR-product 



1(4)1/2 



< 1, 



(IV. 16) 



which reduces to -E pro duct < 0.5, for the choice g = 1. 
In fact, this choice of g is optimal where the inference 
squeezing is very strong. 

We now determine how to optimize the choice of g, 
where the entanglement is weaker. The best choices for 

gJ e D ) and 

Following |35j . 
9 A 2 (JS - gJ° D ) = 2. 9 A 2 J* - (JS, 3%) - (J D , J c ) 



g and g' are adjusted to minimize A 2 (Jq 



dg 



= 



implies the optimal g is given by 

2g=((J c ,J e D ) + (J e D ,J c ))/A 2 J e D . 



(IV. 17) 



(IV.18) 



We note (Jc^Jd) = (Jd'Jc)- There is similarly an opti- 
mum for the value of g', with phase 8' = 9+tt/2. Figure[5] 
shows the optimal value of factor g and g' with different 
cross couplings, versus time, for atoms Na — 200. 

Using the optimal values of the constants, one can eval- 
uate the predictions for the variances, A 2 ( Jq — gJf))- 
We can express the required moments in terms of the 
moments of modes a, 6, e.g., 

Jc ~ 9JE = 9-0 a + Jb) (IV. 19) 

. alb 2 - bla 2 - a\bi + b\ai . 
2 ], 

where we define g± = (l±g) /2, and use 6 2 — (j> 2 = 
0\ — 4>i — due to the symmetry of a and b. The min- 
imum variance is obtained by substitution of optimal g, 
Eq. drVUl. One finds 



where 



A 2 (^ 



2 jZ 



AV. 



9Jd) 



#1 A 2 Jf + 5+AVf , (IV.20) 



b\a 2 



a} Si 



A 2 J$ = A Z (J£ + J£) = A z (Ji + J§). (IV.21) 
Similarly 

(IV.22) 



J, 



x 



< -9Jd =9-(Ja+Jb) 

-g+ 



a\bi — b\ai + a\b 2 — b\a 2 . 




Figure 5. EPR paradox entanglement, shown by the si- 
multaneous inference squeezing of S- = [A 2 (JS- gJ B D )/n ] 

(solid lines) and S+ = \a 2 (J b c +w/2 + g' J* +7r/2 )/n ] (dashed 
lines) with optimal g and g' , where the shot noise level is 
no = (\{Jc)\)/2. Plots show different local cross couplings: 
(a) B = 9.116 G, and (b) B = 9.086 G. Insets show the opti- 
mal value of factors g (solid lines) and g (dashed lines) with 
time. Here Na = Nb = 200, and other parameters are as for 
Fig. [2] 

with 6 2 — (f>i = —{6\ — (f> 2 ) — 7r/2, for which the minimum 
variance is 



A 2 (J*- 5 j£)=. 9 2 A 2 J^+ 5 2 A 2 j5, (IV.23) 



where 



A Z J* = A z (j£ - J£) 

rJi htz. j_ zti ?it 

A 2 ( 



2 . a 2 b\ — b 2 a\ + a[b 2 — b[a 2 , 



Also. 



A 2 J* = A 2 (j£ + J$) = A 2 (jf + J§ ). (IV.24) 



Jc - gJE, Jc - gJ§ > = g- (J+, J x " 



+gl(jf,J*). (IV.25) 
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Figure 6. EPR paradox is predicted (-EEPR-product < 1), for a 
variety of local cross couplings. Here N = Na = Nb = 2000; 
B = 9.116 G (solid line), B = 9.086 G (dotted line), and with 
no cross-couplings, i.e., gi2 = 0, (?22 = gu (dashed line). We 
use optimal value of factors g and g' for each case. Other 
parameters are as for the last figures. 



The minimum interference squeezing A 2 ( J c — gj D ) and 

A 2 (J^, +7r / 2 + g' J^, +7r ^ 2 ) with optimal choices of g and g' 
are shown in Fig. [5j Unlike the original formulations of 
the EPR paradox, in this case, the two inference vari- 
ances are asymmetric. The second inference variance ex- 
ceeds the quantum limit for large enough r. We note 
that for strong correlation as shown by S+ over a large 
range of r the optimal choice for g becomes 1. On the 
other hand, for the poor correlation shown by S— as r 
becomes larger, the optimal choice becomes g' = 0, so 
that the variance S+ is limited to the variance of Bob's 
spin (Fig. [5). 



Figure [6] shows regimes for which the EPR criteria 
are satisfied, i.e., i?EPR.-product < 1- f° r different cross 
couplings. As the number of atoms increases, the spin 
squeezing increases and so too does the degree of EPR 
paradox, as shown in Fig. [7J The result is consistent 
with previous studies of the CV EPR paradox non- 
locality (8j [34] . An EPR paradox for the quadrature 
phase amplitudes of optical modes has been confirmed 
in a number of experiments Whether this effect is 
realizable for atoms however is a different question. Con- 
siderations include the size of r compared to the deco- 
herence time of the BEC condensate and whether inter- 



actions like nonlinear losses ignored in the model (IV.l) 
will come into play to reduce the EPR entanglement. To 
address these questions, we employ the stochastic trun- 
cated Wigner function technique valid for large ./V |33| . 
in the next section. 
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Figure 7. Effect of atom number: (a) N = Na = Nb = 200 
(dashed lines) and N = 2000 (solid lines) on the EPR para- 
dox entanglement with optimal gain factors g and g' , for 
fixed couplings B = 9.116 G. Inset shows the individual 

squeezing inferences S- = ^A 2 (Jq — gJf))/noJ (solid lines), 

S+ = [A 2 (Jc W2 + ff'^ +7r/2 )/ n o] (dashed lines) with opti- 
mal g and g' , for N — 2000. (b) The optimal -EEPR-product 
versus different number of atoms. Inset shows the correspond- 
ing g, g' versus N. 



V. STRATEGY II: SIMULTANEOUS 
EVOLUTION WITH TUNNELING PRESENT 
AND INCLUDING LOSSES 

The degree of EPR entanglement is limited by the 
number of atoms (TV) in the ensemble. Figutre [7] shows a 
base value of ~ 0.83 for TV = 200, but the entanglement 
improves to ~ 0.65 for N = 2000. However, the unitary 
evolution approach of the previous section becomes lim- 
ited once dissipation effects are important. These effects 
are known to occur in ultra-cold atomic systems espe- 
cially at high densities, due to spin-changing atomic colli- 
sions which cause density-dependent two-body and three- 
body losses. With this in mind we investigate the effec- 
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Figure 8. (Color online) Entanglement -E pro duct (a) and Figure 9. (Color online) Entanglement -Eproduct (a) and 
-EEPR -prod uct (b) after evolution of the full Hamilto- ^BPR -produ ct (b) after evolution of the full Hamilto 



nian (III. 2 l including tunneling for N = 200 and B = 9.116 G, 



without the application of the beam splitter. The results are 
shown for values k = 0.1 (blue solid lines), k = 0.8 (red 
dashed lines), k = 0.9 (green dotted lines) and k = 1 (cyan 
dash-dotted lines). 



tiveness of using a large number approximation, namely 
the truncated Wigner function technique |33| . In this 
approach, higher order terms in 1/N 3 / 2 are ignored, to 
allow a stochastic calculation based on a positive Wigner 
function. This method readily scales to large numbers of 
atoms and modes |49| . and can include nonlinear losses. 
The detailed description of the method can be found in 
Appendix ^] 

The predictions of the truncated Wigner method are 
indistinguishable from the exact method given in Ap- 
pendix [A] in the case of zero losses and no tunneling, 
which confirms its validity for N - 200-2000. The 
advantage of the Wigner method is that it allows a 
ready solution of the dynamics of the full Hamilto- 
nian ( III. 1 III. 2) where both tunneling and nonlinear 
terms are present, even for large atom numbers. The 
method also allows the inclusion of losses, which are 



nian (III. 2 1 including tunneling for TV = 2000 and B = 
9.116 G. Meaning of the lines is the same as in Fig. [8] 



known to destroy squeezing and entanglement, and will 
come into play in realistic experimental arrangements. 
For example, tunneling cannot be completely suppressed, 
and also potentially significant is the role of nonlinear 
loss, that will come into play with large atom numbers. 

Once tunneling is present in the nonlinear Hamil- 
tonian, entanglement can be created between the two 
modes, without the second beam splitter step described 
in Section |IV| Figure [8] shows the effect of tunneling on 
the entanglement. Strong tunneling (k ~ 1, i.e., with the 
strength of the same order as the nonlinear interaction) 
produces significant entanglement even without the final 
application of the beam splitter. However, tunneling this 
strong is hard to achieve in a simple two- well BEC exper- 
iment, where values of k ~ 10 -3 . . . 10 -2 are more com- 
mon. The achievement of such large couplings is likely 
to require a more sophisticated experimental design. 

If one increases the number of atoms N, with the 
absence of losses the dimensionless drift part of equa- 
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Figure 10. (Color online) Ent anglem ent -E pro duct after evo- 
lution of the full Hamiltonian (III. 2 1 including tunneling for 
N = 200 and B = 9.116 G, (a}~without and (b) with the 
application of the beam splitter. The results are shown for 
values k = 0.01 (blue solid lines), k = 0.5 (red dashed lines), 
and k = 1 (green dotted lines). 



tions (B.4) will stay the same, but the results will differ, 
as shown in Fig. [9] The entanglement improves with 
higher N, as in the case of the two-step strategy and 
as found previously for ground state calculations. This 
result is consistent for continuous variable EPR entan- 
glement, which has been predicted in optics for high in- 
tensity Gaussian states. 

We also analyze the predictions for entanglement with 
the insertion of step 2, the "beam splitter" interaction 
described in Section |TV] B, the objective being to under- 
stand whether a small amount of tunneling and nonlinear 
loss present in the first interaction stage will detract from 
the amount of entanglement that is predicted by the two- 
step strategy. We focus on the case where local couplings 
are present, although this gives a worse case prediction 
for entanglement. The losses and tunneling tend to de- 
tract from the amount of entanglement, so that in the 
case of 2000 atoms, with local couplings, an entangle- 



Figure 11. (Color online) EPR paradox entanglement as mea- 
sured by -EEPR-product without tunneling and with the appli- 
cation of the beam splitter, for B = 9.116 G and Na = Nb = 
2000. (a) Only linear losses are enabled; no losses (solid blue 
lines), 71 = 5 x 10 -3 (dashed red lines), 71 = 10~ 2 (dotted 
green lines), (b) Only two-body intra-species losses are en- 
abled; no losses (solid blue lines), 722 = 5 x 10~ 6 (dashed red 
lines), 722 = 10~ 5 (dotted green lines). 



ment of ~ 0.4 (enough to give an EPR paradox) as shown 
in Fig. [4] is reduced to 



0.6 (Fig. 10 1. We emphasize 
here that losses are highly adjustable through changes 
in density, so these issues are more related to appropri- 
ate experimental design than to fundamental limits. The 
application of the beam splitter improves entanglement 
significantly when tunneling is weak, but with k values 
close to 1 tunneling becomes the prevalent source of en- 
tanglement, as illustrated in Fig. 
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The presence of linear losses, or intra-species losses de- 
creases the maximum entanglement, as shown in Fig. |11| 
But when the nonlinear inter-species losses are enabled, 
our simulations show that they unexpectedly increase the 
entanglement (Fig. |12[). 

We have found that this only occurs with a nonlinear 
loss 712 that specifically couples the two species together. 
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Figure 12. (Color online) EPR paradox entanglement as mea- 
sured by -EEPR-product without tunneling and with the ap- 
plication of the beam splitter, for B = 9.116 G and Na = 
Nb = 2000. Only two-body inter-species losses are enabled; 
no losses (solid blue lines), 712 = 10~ 4 (dashed red lines), 
712 = 10~ 3 (dotted green lines). 



loss, we calculate that enhanced entanglement is possible. 
Finally, we note that the optimum regime is for rather 
large coupling or tunneling values between the spatial 
modes, which appears to require a different experimental 
design to a simple two-well system. 

The EPR paradox entanglement studied in this paper 
confirms an inconsistency of local realism with the com- 
pleteness of quanutm mechanics |7J, and is evidence for 
the form of nonlocality called "steering" [lOj . We make 
the point however that EPR entanglement is not itself 
sufficient to imply a direct failure of local realism, as 
would be demonstrated by a violation of a Bell inequal- 
ity (6j [8] . The method employed to arrive at the predic- 
tions of EPR entanglement illustrates this point, since 
the truncated Wigner function is positive, and therefore 
provides a local hidden variable theory to describe the 
statistics of experiments where the measurement is effec- 
tively that of a quadrature phase amplitude, even though 
significant EPR entanglement can be obtained [8]. The 
distinction between the EPR steering and Bell forms of 
nonlocality has been emphasised further in recent papers 



Other forms of loss, including linear loss, will simply re- 
duce the entanglement — as one expects from decoher- 
ence effects. Generating entanglement by a manipulation 
of reservoirs has been observed for atomic ensembles |50j . 
The process in our case appears to be related to the fact 
that nonlinear absorption just by itself is known to create 
a nonclassical state |51j . which can then become entan- 
gled through linear couplings alone. 
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VI. SUMMARY 



Appendix A: Exact dynamical solutions 



In summary, we have analyzed two strategies for dy- 
namical preparation of squeezed and entangled atomic 
states, through the use of nonlinear interactions that rely 
on the naturally occurring S-wave interactions between 
trapped atoms. The possible advantage of the meth- 
ods proposed here is that they do not require the use 
of separate local oscillators, as often employed in optics 
for measuring entanglement. These could be potentially 
technically difficult to use in some cold-atom systems, 
due to dephasing that is caused by interatomic interac- 
tions in the local-oscillator itself, when combined with 
number fluctuations. However, atomic homodyning has 
recently been realised [3 , and may present a better strat- 
egy in some cases, particularly where it is important to 
avoid local cross-couplings. Work done in parallel with 
ours studies the dynamical prepration of entanglement 
for this case, using a single mode at each well |52| . 

We find that robust spatial entanglement and EPR 
inference is possible under the correct conditions, even 
including Poissonian number fluctuations. These effects 
are also not greatly perturbed by realistic loss values. In 
fact, with the presence of certain types of inter-species 



a. One-mode model: 

We start by considering a simple nonlinear, single- 
mode interaction that occurs locally at each well, and 
is modeled by Hamiltonian 



H/h= y -: N 2 



(A.I) 



In order to understand the dynamics induced by the 
above Hamiltonian, we calculate results in the Heisenberg 
picture, where one obtains: 



lit 



H,CLi 



iga/aa — —igNa. (A. 2) 



Since the number of particles is conserved, this has the 
solution: 



a (t) = exp —igNt a (0) , 



(A.3) 



12 



which gives 

<a(i)) = <¥(0)|o(t)|¥(Q)) 



oo 



C* n C m (n| exp 



-igNt 



o(0) |m) 



second, and that the condensate mode is give by: 

00 

l*(o)) = |o) 01 H 02 = |o} oi ^c n H a2 , (A. 

n=0 

which gives 



00 

E 

n=0 



C* C„+i y/n + 1 exp [— ip/it] 



a| 2 fe" l9t - 



2r> 



afci |tf (0)) = 0, 
4a 2 |*(0)) = H 2 |*(0)) 



(A.9) 



a 



exp [I 



1)] 



(A.4) 



In this coherent state the number fluctuation variance is 
N. Suppose we apply a linear beamsplitter, then, after 
the beam-splitter: 



This predicts a well-known behavior with three char- 
acteristic time scales: 

• On very short time scales, there is simply an oscil- 
lation with a renormalized frequency of u' = gn. 



1 



a 1 



(61 + a 2 ), 



a-2 = y|(" 2 - 



where n 
sically. 



I a I 2 . A similar result is obtained clas- the state becomes 



• On intermediate time-scales, there is a quadratic 
damping, with a characteristic damping rate of 
gV^, which is proportional to the standard devi- 
ation in the initial particle number. 

• Finally, on very long time-scales there is a suc- 
cession of periodic revivals where the initial state 
and all its properties are regained exactly, apart 
from a possible phase-shift. This occurs whenever 
t = 2k /g . 



b. Two-mode model: 

Next, consider a simple two-mode model, with two in- 
ternal (spin) modes, where the Hamiltonian for the cou- 
pled system is: 

& = 2 E^H^ = 2 '-^SijNiNj : . (A.5) 

We can solve this using either Schroedinger or Heisen- 
berg equations of motion. In the Heisenberg case, one 
obtains: 



ddi i 
~dt ~ h 



H,di 



(A.6) 



Since the number of particles is conserved in each mode, 
this has the solution: 



hi (t) = exp 



a % (0) . (A.7) 



We suppose the initial quantum state factorizes into 
a vacuum state in one mode and coherent state in the 



l*(0)> = 



V2 



(A. 10) 



(A.ll) 



After the application of the nonlinear Hamiltonian 
(A.5 1, squeezing is generated locally, and 



(oi(t)) = (*(o)|a,(t)|*(o)) 



V2 V2 



-i^gijNjt 



V2' 



E 

ni ,ri2 ,mi ,7712=0 

x(ni\(n 2 \ exp 



V2 



exp 



2 



-i^gijNjt 
3 

( e -igut _ ^ 



ai\mi)\m 2 ) 



exp 



(e- iBat - 1) 



(A.12) 



We wish to calculate the phase variance in 9 = 6\ — 6 2 , 
as this gives rise to the decay in an interference pattern: 

((A§) 2 ) = (8 2 ) - (§) 2 . (A.13) 



Appendix B: Truncated Wigner Method 

Including losses, the master equation in four-mode ap- 
proximation is written as 



dp 



l -[H 7 p]+Y,^i[pl (B.l) 



where the Hamiltonian is defined as (III. 2 ), and loss term 
has the form 



Ci\p] = 20 lP 6\ - OjOip - pOJOi. (B.2) 
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We consider three different sources of losses: O22A = a 2 
and O22B = b 2 , (two-body intra-species loss), 0\2A = 
a\a,2 and 0\2B = (two-body inter-species loss), 

0\a — a,2 and 0\b — b\ (linear loss). 



In order to further normalize the equation ( B . 1 1 , we 
use the dimcnsionlcss time r = guN^t introduced in 
Section |ITT} This gives us the dimensionless master equa- 
tion 



dp 
Jr 



i[P\ 



(B.3) 



with a dimensionless Hamiltonian H as in Eq. (III.3I, 
and dimensionless loss coefficients 7j = Ji/quNa- This 
equation can be transformed to the equivalent partial 
differential equation by applying the Wigner transforma- 
tion 1351: 



IV A = — I d 2 X 1 d 2 X 2 d 2 X 3 d 2 X 4 

Y[exp(-\iZ* +X*z l ) 
\i=i ) 

x Tr I A n ex p ( a 4- k^) I , 

where 4-vectors z T — (a± Pi a.2 P2) and z T 
in bi 0,2 62) were introduced for convenience. The 
resulting differential equation, after truncating higher- 
order derivatives, is a Fokker-Planck equation (FPE) for 
the truncated (positive) Wigner function W = W [p] 



x 



dW 
~dV 

where 



= -d T z aW - dJ,a*W + Tr {d z d T z ,BB H } W, 



z = (at Pi a 2 P2) , 



/K1P1 + ai{gn|ai| 2 +3i 2 |a 2 | 2 }N 
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This FPE is equivalent to the following set of stochastic 
differential equations (SDEs) |53| : 



dz = adr + BdZ, 



(B.4) 



where dZ T = {dZ\2A dZi2B ■ ■ ■ dZ\A dZ\B) is a com- 
plex 6-dimensional Wiener process. These equations can 
be solved numerically using conventional methods, and 
their solution can, in turn, be used to get the expecta- 
tions of symmetrically ordered operator products as 



sym 



<«r<. 



..Wd 2 z 

) paths 5 



where () pa ths stands for the average over the simulation 
paths. 

As for Section [m] we assume the initial state to be the 
coherent state: 



W) = \ao)a 1 \ao)a 2 \Po)b 1 \Po)b 2 , 
where a = ^N A /2, /3 = y/N B /2. Therefore 



(B.5) 



qT I _d_ Odd 

u z — \ 9ai 9/3i da 2 3/3 s 



-tadrift — aj 



Oii = a + -7]u, Pi = po + -r) 2i , 

where r/u and r)2i are complex normally distributed ran- 
dom numbers with (77^77*./) = SjkSu. 
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